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Abstract. We derive, and experimentally demonstrate, an intcrferometric scheme 
for unambiguous phase estimation with precision scaling at the Heisenberg limit that 
does not require adaptive measurements. That is, with no prior knowledge of the 
phase, we can obtain an estimate of the phase with a standard deviation that is 
only a small constant factor larger than the minimum physically allowed value. Our 
scheme resolves the phase ambiguity that exists when multiple passes through a phase 
shift, or NOON states, are used to obtain improved phase resolution. Like a recently 
introduced adaptive technique [Higgins et al 2007 Nature 450 393], our experiment 
uses multiple applications of the phase shift on single photons. By not requiring 
adaptive measurements, but rather using a predetermined measurement sequence, the 
present scheme is both conceptually simpler and significantly easier to implement. 
Additionally, we demonstrate a simplified adaptive scheme that also surpasses the 
standard quantum limit for single passes. 
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1. Introduction 

The interferometric measurement of an optical pliase is an indispensable tool in 
metrology, used in length measurement, optical characterization and a myriad of other 
applications. Phase measurement is closely related to the more general concept of 
quantum parameter estimation, which is important for characterizing Hamiltonians 
for applications in quantum information technology, for example. In practice, there 
exist many different types of phase measurement problems, each with its own set 
of requirements and limitations, which has led to a diverse range of interferometric 
techniques. 

In this paper, we consider the estimation of a completely unknown phase G [0, 27r) 
in the situation where the total number of photon-passes through the system being 
probed is the crucial resource, and there is no inherent limitation to the number of 
times each photon can pass through the phase shift [H [21 |3l H] (see Appendix). This 
corresponds to the case of a single measurement of a fixed phase. The number of photon- 
passes is important because, for example, the light may cause damage to a sensitive 
sample, or each sampling of the phase may be difficult to achieve (such as in probing 
the strength of a Hamiltonian coupling for a quantum information application [5]). In 
this regime, the aim is to obtain an estimate (^est with minimal uncertainty A0est for 
a given number of applications, N, of the phase shift [31 H]. In the nomenclature of 
quantum information, this restriction equates to having a limited number of queries of 
the system. 

The regime we consider may be contrasted with the "phase sensing" regime, such as 
would be found in the context of gravity-wave detection or atomic clocks. In that case 
the aim is to measure changes in a time-varying phase. There is an inherent limitation 
to the number of passes of each photon, because additional passes will eventually lead 
to less accurate results due to the varying phase. In addition, because the phase shifts 
are small, the phase is already approximately known, allowing any phase ambiguities 
to be resolved. In that regime, particular techniques and states of light, such as n- 
photon path-entangled "NOON" states [6l[7], are desirable due to their improved phase 
resolution from a single pass through the sample. The phase ambiguity in the NOON 
state (which only allows the phase to be estimated modulo 27r/n) is resolved using the 
pre-existing approximate knowledge of the phase. 

In phase measurement, there are two well-known limits to measurement accuracy. 
One is the standard quantum limit (SQhj^l, which is the limit obtained by "standard" 
techniques, such as coherent states and photon counters. The other is the Heisenherg 
limit, which is the fundamental limit using arbitrary quantum states and measurements. 
When there is an inherent limitation to the number of passes that may be used, the SQL 
and Heisenberg limit are 0(1/ a/w) and 0{l/n), respectively, where n is the number of 
photons [HI |9]. To attain the Heisenberg limit in this case it is necessary to create 

X In some areas this is called the shot- noise limit, with the term "standard quantum limit" reserved 
for a limit set by additional constraints which are not relevent to the situation we consider. 
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path-entangled n-photon states, such as NOON states [HI Ej or other non-classical 
states [HI, [9l [TOl [11], [121 HSl El- However, such states are very difficult to produce 
experimentally even for moderate n[T5lfT6l[T71[T8[[T9l[2n[[2Tt[22]. 

In this paper, we follow some recent papers [Tl[2] in using the term Heisenberg limit, 
as defined above, for the general situation, regardless of whether there is a limitation on 
the number of passes of any individual photon. In terms of the total number of photon- 
passes, N, the Heisenberg limit is 0(1/A^) [2, [3]. In order to give a meaningful basis 
for comparison, we take the SQL to be the maximum accuracy possible with classical 
states and single passes, which is 0(l/v^)- When there is no inherent limitation to the 
number of passes for each photon, interference fringes with a width scaling as 0(1/A^) 
may be obtained simply by using N passes of a single photon. The difficulty is that 
this also yields an ambiguous phase estimate, in exactly the same way as using NOON 
states. When there is no additional knowledge of the phase, resolving this ambiguity is 
a major problem. 

Recently, a number of approaches have been proposed to resolve the above 
ambiguity, both in the context of phase estimation [U [3] and in related contexts [21 
[231 [21|- These are closely related to the quantum phase estimation algorithm (QPEA) 
of Cleve et al (221 [26] (which is the heart of Shor's factorization algorithm [2Z]), as well as 
Kitaev's phase estimation algorithm [21]. None of these protocols were rigorously shown, 
either numerically or analytically, to give Heisenberg-limited scaling for the standard 
deviation of the phase estimate (see section [2]). Another scheme, not related to the 
QPEA, was shown numerically to be able to achieve Heisenberg-limited operation [29j . 

Heisenberg-limited phase estimation was realized experimentally for the ffist time 
in [30] , using a QPEA-inspired adaptive phase estimation scheme proposed in that work. 
That this protocol does achieve Heisenberg-limited scaling was also shown using rigorous 
numerics. In fact, two schemes were considered in [30]. Firstly, the QPEA itself was 
implemented, showing that this scheme fails to achieve a phase uncertainty below \h/N 
due to the high wings on the distribution. Secondly, a generalization of the QPEA[§[ was 
presented and demonstrated to achieve an uncertainty within a factor of 1.6 of the 
Heisenberg limit, which approaches ti /N for large N [TTl [131 [311 [32] ■ Subsequently, a 
different QPEA-inspired adaptive phase estimation scheme was analytically shown to 
give Heisenberg-limited scaling by Boixo and Somma [1]. 

Higgins et al [30] demonstrated the QPEA and its generalization using multiple 
passes of a single photon through the phase shift. However they also noted that the same 
schemes could instead be implemented equivalently using single passes of NOON states. 
An advantage of using NOON states is that the time taken for a measurement does not 
scale with N , thus allowing the Heisenberg limit to be attained in the phase sensing 
regime. A disadvantage of using NOON states, or related multiphoton entangled states, 
is that the size of the state scales as N . Another disadvantage is that the NOON- 
state protocols require photon- resolving detectors with inefficiency at most 0(1/A^). By 

§ Note that in [3D], the (generahzed) QPEA was called the (generalized) Kitaev algorithm. 
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contrast the multipass protocols can achieve Heisenberg-limited scaling for any non-zero 
detection efficiency. See the end of section [5] for further discussion. 

Here we present and experimentally demonstrate two new phase estimation schemes 
that offer significant simplifications compared with the generalized QPEA of [30], while 
still giving accuracy better than ©(I/a/ZV). Both can be applied to either NOON states 
or multiple passes of single photons, as with the QPEA and its generalization. As 
in [30], we have performed the experimental demonstration using multiple passes of 
single photons, as shown in figure [TJ Whereas the Heisenberg-limited scaling of the 
generalized QPEA of [30] was shown numerically, we here provide analytical proofs 
bounding the scaling of our schemes. 

The ffist scheme, which we shall call the QPEA-standard hybrid or simply 
hybrid scheme, combines the original QPEA with standard (single-photon, single-pass) 
interferometry to achieve a phase uncertainty 0((ln A^)^/^A^~^/^) or better. The second, 
non-adaptive scheme uses a predetermined sequence of single-photon measurements. 
Remarkably, this non-adaptive scheme achieves a phase uncertainty 0(l/iV), although 
with a higher overhead than the adaptive scheme of [30j. As well as being of fundamental 
interest, the reduction in complexity offered by these new schemes makes the practical 
adoption of multipass phase estimation techniques more attractive. 

2. Existing schemes 

In order to compare the performance of phase estimation schemes, it is necessary 
to define an appropriate measure of uncertainty. Uncertainty measures are not 
interchangeable. For example it may be the case that although an accurate phase 
estimate is obtained the majority of the time, a small proportion of measurement errors 
are large [I3|. This will lead to a large standard deviation in the phase estimate, 
although other measures (e.g. a confidence interval) may fail to capture this uncertainty. 
We use the standard deviation to unambiguously quantify phase uncertainty — if the 
standard deviation of the phase estimate scales as 0(l/iV), so will other measures of 
uncertainty [33]. The converse is not generally true. 

When considering a single measurement of an unknown phase, it is necessary for 
the measurement procedure to be self-contained in order to fairly evaluate it. That 
is, the procedure should not use additional phase information from some other source 
without counting the resources used to obtain that information. There are two main 
issues to be considered here: 

(i) NOON states or multiple passes through a phase shift provide an ambiguous phase 
estimate, and extra measurement steps are required to resolve the ambiguity. 

(ii) Additional information about the phase may also be useful during the measurement 
procedure to design an optimal or near-optimal measurement, as is the case for 
measurements based on squeezed states, for example [TO] . 

One approach to obtain a self-contained, unambiguous measurement is to use the 
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Figure 1. Schematic of the experiment. A polarizing beam splitter (PBS) creates 
high-fidelity equal superpositions of left- and right-circular polarization modes, which 
form the two arms of the interferometer. The 2''9 half- wave plate (HWP) in a motorized 
rotation stage implements the controlled phase shift. A photon passes multiple times 
through the HWP, implementing the multiple applications of the unknown phase 
shift. Quarter-wave plates (QWPs) undo 7r-phase shifts in the horizontal/ vertical 
polarization basis upon reflection. Mirrors on motorized translation stages select 
the number of passes the photon undergoes. The photon is then discriminated 
by a polarizing beam displacer (PBD), filtered to remove unwanted background, 
and detected by a single-photon counting module (SPCM). In the hybrid scheme, 
measurement results determine the setting of the controlled phase shift as indicated 
by the dotted path from the processor to the control phase box. This path is not used 
for the non- adaptive scheme. 



QPEA [251 126]. The QPEA involves two steps. In the encoding step, K + 1 qubits, each 
originally in state (|0) + \ l))/\/2, control the application of a sequence of phase shifts 
exp(i2^0), exp(i2^~-'^0), . . . ,exp(i0). This prepares the state Z]n=o s^"*^!^), where \n) is 
a state in the logical basis, and N = 2^^+^ — 1. In the measurement step, the inverse 
quantum Fourier transform is used to perform a measurement in an entangled basis on 
the encoded state. This measurement is the optimal one for obtaining an estimate 0est 
of the unknown phase (j) encoded in the state p6l [25]. Remarkably, the measurement 
step can be implemented without entangling gates by using a Markovian sequence of 
local adaptive measurements on each of the K+1 qubits. At each step the measurement 
basis depends only on the previous measurement and its result [3^ . 

Although the measurement step is optimal for phase estimation, the encoding step 
is not. As a consequence, the phase uncertainty A0est = (where Vh is the Holevo 
variance [32j) is dominated by the high wings of the distribution of phase estimates, and 
the QPEA does not achieve a phase estimate with uncertainty scaling at the Heisenberg 
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limit 71 /N. Rather, the uncertainty actually scales as I/a/ZV [3U]. 

To achieve scaling at the Heisenberg limit, the QPEA can be modified by repeating 
each application of the phase shift operator exp(i2'^0), k G {0, . . . , K}, on each of 
M independently prepared qubits [30]. Bayesian adaptation (that is, controlhng the 
measurement made on each qubit based on the distribution for conditioned on all 
preceding measurements [13]) can replace the single qubit measurements to achieve 
0(1/A^) scaling for M > 4. 

As explained above, a phase shift of p(f) can be achieved by using an interferometer 
and either a single photon with p passes through the unknown phase shift, or a NOON 
state with n = p photons and a single pass. In the case of multiple passes, the probability 
distribution for the phase based on the measurement is P{(f)\±) = (l±sin(p</)))/27r, where 
the sign depends on which output the photon is detected in. In the case of NOON 
states, all information about is again contained in a binary outcome ±, this time 
corresponding to the parity of the photon number at one of the outputs. Moreover, 
the probability distribution for the phase is identical to that obtained in the single- 
photon case. It is evident that both cases only yield information about (p modulo 27r/p; 
as explained, the role of the QPEA, or the generalized QPEA scheme in [30], is to 
resolve this phase ambiguity. The experimental demonstration in [30] was achieved 
with multiple passes. 

The total resources used is quantified by the number of applications of the phase 
shift, and is = M(2^+^ — 1). When using multiple passes of a photon this is equal 
to the total number of photon-passes, and it is equal to the total photon number for 
the equivalent scheme using NOON states. As discussed in [3] and the Appendix, this 
is the appropriate way to quantify resources for problems of this type. 

3. QPEA— standard hybrid scheme 

For M = 1 the Bayesian adaptive scheme of [30] reduces to the QPEA. However, as 
noted above, the control protocol in the M = 1 case has a much simpler formulation 
as a Markovian scheme. If the problem of outliers of (pest (due to the high wings of 
the distribution) could be solved, the QPEA would be an attractive scheme due to its 
simplicity. 

Here we present a solution, resulting in an uncertainty that beats the SQL, in which 
the QPEA is supplemented by some additional simple measurements. These additional 
measurements are those of standard interferometry, comprising single phase shift 
applications on many single qubits without adaptive measurement. A related scheme 
was explored in [29], using standard interferometry to supplement adaptive NOON state 
inputs with adaptive measurements; numerical simulations showed Heisenberg-limited 
scaling with an overhead factor of 2.03 [29]. We now analytically determine the scaling 
for the phase uncertainty using the adaptive scheme proposed above. 

Consider using resources to implement a standard phase measurement, in 
addition to the QPEA using A^^q = 2^^^ — 1 resources. First we consider the uncertainty 
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in the standard phase measurement, and for simphcity alternate the controllable phase 
shift 6 (see figure [1]) between = and 62 = tt/2. There are thus two sets of 
measurement results governed by probabilities Pi = (l + cos0)/2 and P2 = (l + sin0)/2 
respectively. For even iVg, we obtain estimates of Pi and P2, denoted z/i and 1^2, each 
based on Ns/2 measurements. The estimate of is then determined from the phase of 
{2i'i — 1) + i(2z/2 — 1). Using the Chernoff bound |35] (or Hoeffding's inequality [36]). 
the probability that the errors in these estimates exceed e are bounded as 

P{Wk -Pk\>e)<2 exp {-2Nse') , (1) 

for k G {1,2}. The probability that either of the estimates differs from the P^ by at 
least e is then upper-bounded by 4exp (— 2iVse^). Using simple geometry, it is easily 
seen that ii {uk — Pk\ < ^, then the error in the angle can not exceed arcsin {2y/2e). Let 
us take 



/3/(iVs)/32iVs, (2) 
where / satisfies f{x) < x for x > 1. With this choice, we find that 



arcsin = arcsin J3f{Ns)/4:Ns 



< in/3)^f{Ns)/Ns. (3) 

In the second line we have used the fact that Ns > 1, so 2\/2e < v^/2, and the convexity 
of the arcsin function over the range [0, -1/3/4]. Therefore, the probability of the phase 



error being at least 6(f) = {n/S)^ f{Ns)/Ns is no greater than 4exp [—3f{Ns)/lQ]. 

We now wish to place bounds on the phase uncertainty. We quantify the phase 
uncertainty by the square root of the Holevo variance, Vn [37j, where 

V^H = |(e^('^"^^='))P^ - 1. (4) 

When small, the Holevo variance is well approximated by [13] 

l^ = 4(sin2[(0-0est)/2]). (5) 

To place an upper bound on this variance, we consider a suboptimal method of 
determining the phase estimate. Consider combining the standard phase estimate, 0s 
and the QPEA phase estimate, 0q, by using (pq when the estimates differ by less than 
2S(f), and using 0s if the two estimates differ by more than 250. There are then three 
different alternatives which need to be considered. 

(i) If the standard phase estimate is inaccurate, in the sense that |0s — 0| > 50, then 
the contribution to V can not exceed 16exp [— 3/(Ai's)/16]. Choosing f{Ns) = 
min[A^S) (32/3) In Ai's] ensures that this contribution is 0(A^g'^). 

(ii) For the case where the standard phase estimate is accurate (|0s — 0| < 50), if the 
QPEA estimate differs from the standard estimate (so |0s — 0q| > 250), then the 
QPEA estimate must be inaccurate (so |0q — 0| > 50). Provided the QPEA uses 



Heisenherg -limited phase estimation without adaptive measts 



8 



a randomized initial setting of the controllable phase, the probability distribution 
is 

PU ^ - sin^ [(iVQ + l)(0Q-0)/2] 
'^^'^^^ 2vr(iVQ + l)sin^[(0Q-0)/2]- 
Using this probability distribution, the probability of |0q — 0| being greater than 



(50 does not exceed [9/{Tr'^NQ)]^ Ns/ f{Ns), and the overall contribution to V then 
does not exceed 

(l/iVQ)v7(iVsj7iVs. (7) 

(iii) Next consider the case that the standard phase estimate is again accurate 
(10s ~ 01 < '50) and the standard and QPEA estimates agree (so |0s — 0q| < 
2(50). This means that the QPEA estimate must be accurate, in the sense that 
|0Q — 01 < 3(50. With this limit, we find that the possible contribution to V is 
bounded by 

i2/NQ)^fiNs)/Ns + 2/inNl). (8) 
Collecting together all three cases gives the upper bound on V as 




+0iN,') + 0iN^'). (9) 

Ignoring the slow variation of f{Ns), the minimum of ([9]) as a function of = A'q + A^g 
will be obtained for A"s = A^/3 and A'q = 2N/3, which gives 

V = 0{Vh[NN-^/^). (10) 

We have implemented the above hybrid scheme both experimentally (see below) 
and numerically. We do not consider special cases as in the above derivation, but 
rather determine the best estimate from the normalized product of the distributions 
conditioned by Bayesian updating. Also, we increment the controllable phase 6 by 
7r/A"s after each measurement, rather than alternating between and tt/2. Numerically 
we find that choosing A"s = 2^ , so that A^ = 3 x 2^^ — 1, is close to the optimal choice, 
as expected. Moreover, we find that for this choice Vh x N^/'^ is nearly constant, also 
as expected (it increases from 4.83 at A^ = 5 to 6.17 at A^ = 767). 

It is useful to contrast our hybrid scheme with the generalized QPEA of [30]. In 
the latter, M (the number of independent iterations for each 2'^-fold application of 
the phase shift), did not depend on k. In our hybrid scheme M depends on k, as 
M{K, k) = 1 + 5kfl2^ . The hybrid scheme also differs in using mostly non-adaptive 
measurements {2^ + 1 non-adaptive measurements versus K adaptive ones). These 
observations suggest that, by considering other dependencies M{K,k), good schemes 
may be found that dispense with adaptivity altogether. 

4. Non-adaptive scheme 

In the adaptive case, it is possible to obtain scaling at the Heisenberg limit for 
M > 4 [30j. In contrast, for non-adaptive measurements, using small values of M 
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such as these results in much larger variances that do not decrease below some fixed 
amount (see figure [2]). Rather than there being a threshold value of M above which 
Heisenberg scaling is obtained, these results suggest that M must increase as a function 
of K to avoid results with large error. But this can not yield Heisenberg scaling because 
it leads to an increasingly large overhead. 



M = 4: (numeric calculation) ^^^^ 

M = 6 (numeric calculation) 

M = 8 (numeric calculation) 

-- ■- M = 12 (numeric calculation) ^^^^ 

Standard quantum limit (numeric calculation) 

Heisenberg limit (analytic calculation) ^^^^ 

I I I I I I ' ... I'M 

10 100 1000 lO'* 

Number of resources, A'^ 

Figure 2. Predictions of standard deviations A(/)est of distributions of phase estimates 
under a non-adaptive scheme with fixed M and controllable phase 6 incremented by 
tt/M after each measurement. Numerical predictions were based on 2^** samples. The 
Heisenberg limit and the SQL are shown for comparison. 

On the other hand, as we will now show, with a suitable choice of M{K, k) it is 
indeed possible to obtain scaling at the Heisenberg limit without adaptive measurements. 
Specifically, we choose M{K, k) = Mk + — k), where Mk is the value of M used 
for the 2'^-fold application of the phase shift, and /i > 0. The optimal performance 
is obtained by numerically optimizing both of these parameters, which we consider 
later. First, we prove rigorously that scaling at the Heisenberg limit is guaranteed for a 
particular choice of Mk and /i. 

Since our scheme is non-adaptive, it can be performed in any order; we consider 
using M{K, k) independently prepared qubits undergoing the 2^-fold application of the 
phase shift in the sequence = 0, 1, 2, . . . , K. Consider the results for A; = 0; applying 
the analysis using the Chernoff bound as in the hybrid scheme, for even M{K, 0) choosing 
e = ^3/2/4 localizes the true phase to an open segment of angular size 27r/3 with 
probability at least 1 - 4exp [-3M{K, 0)/16]. 

Next, using this analysis on the measurements with a two-fold phase shift (A; = 1), 
we see that the phase is locahzed to two diametrically opposite segments of size vr/S. As 
the separation of these segments is 27r/3, and they are open, only one of them overlaps 
with the segment obtained for k = 0. Hence, combining the results of these two sets 
of measurements localizes the system phase to within a segment of size tt/Q with high 
probability. This analysis is continued up to the 2'^-fold phase shift. At stage k the 
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phase is localized to a segment of size (27r/3)/2'^, with probability of error bounded 
above by 4exp [-3M{K, k)/16]. 

Consider the maximum error for each case where the phase may lie outside the 
segment. The maximum error in the case for = is vr, and occurs with probability no 
greater than 4exp [— 3M(i^', 0)/16]. The contribution to ^sin^ [{4>est — 0)/2]^ therefore 
does not exceed 4exp [—3M{K, 0)/16]. Then, for k = 1, . . . , K, the error can not exceed 
{Att/3)/2^ (for example, if the phase is outside the k = 1 segment, the error can be as 
large as the size of the k = segment). The probability of this error does not exceed 
4exp [—3M{K, k)/lQ]. Thus the contribution to ^sin^ [(0est ~ 0)/2]) from these errors 
does not exceed 

4^-3M(X,0)/16^4^3i^2 
k=l 

Finally, if the phase is correctly localized in the k = K segment, the maximum 
error is (27r/3)/2^. Including this contribution we obtain the upper bound on V as 

V<V^^ = iQe-^MiK,o)/w _^ (2^/3)22-2X 

+ 16(27r/3)2 J2 2-2V3^^(^''=)/i^ (12) 

k=l 

The total number of phase shift applications is = Y.k=o M(K, k)2'^. It is evident 
that, for M{K,k) independent of k and K, V^ax = 0(1), in agreement with the 
numerical results. If we allow M to be a function of K but not of k, then the best 
scahng of Vmax is obtained for M oc K, which yields \4iax = 0(ln(A^)/A^^). To determine 
a variation of M{K, k) that does yield scaling at the Heisenberg limit, we may take 
5Af(i^,fc)(K.axA^') = 0, which gives exp[-3A;ln2-3M(ir,fc)/16] = {3/7r^)V^,jN . To 
satisfy this, M{K, k) should have the linear variation M{K, k) = Mk + 16 (In 2){K — k). 
With this choice, it is easy to show that 

Knax < (27r/3)'(l + 32e-=^^^/i6)2-2^. (13) 

Since = 0(2-^), this proves that the phase uncertainty scales at the Heisenberg 
limit when we use linear variation of M{K,k) with fi = 16 In 2. However, A0est has a 
substantial overhead (about 54 for Mk = 23). 

The above result represents an upper bound on the phase variance. However, 
using numerical methods we find that the optimal values of n and Mk are smaller, and 
yield much smaller overhead. We performed numerical simulations (up to = 10^), 
incrementing the controllable phase 6 by TT/M{K,k) after each measurement, instead 
of alternating between and tt/2. These indicate that values of Mk = 2 and fi = 3 
consistently and robustly give a low overhead, less than 2.03 times the Heisenberg limit. 
This is only marginally above the overhead of 1.56 for the adaptive scheme (M = 6) 
of [30], and very close to the scaling constant for the adaptive scheme of [29] . 



(27r/3)/2^ 



-3M{K,k)/lG 
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5. Experiment 



We demonstrate both the hybrid and non-adaptive schemes introduced above using 
photonic qubits in a common-spatial-mode polarization interferometer, as in figure [H 
similar to that of [30]. As with the QPEA and its generalization, these new schemes 
may be equivalently performed either using multiple passes of single photons or single 
passes of NOON states — here we use the former. The two arms of the interferometer 
are the right- and left-circular polarization modes, with phase shifts applied using half- 
wave plates. Computer-controlled motorized stages adjust the configuration of passes 
through the system phase and the setting of the control phase in the other arm of the 
interferometer. Pairs of 820 nm single photons are generated by a type-I spontaneous 
parametric down- conversion source, pumped by a continuous- wave laser diode. One 
photon of the pair is detected immediately, the other is sent through the experiment, 
and a coincidence detection heralds a successful measurement. 
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Figure 3. Performances of phase estimation schemes with varying numbers of 
resources, N. In contrast to figure [51 here we show standard deviations A^cst of 
distributions of phase estimates scaled by ; in this case the SQL asymptotes 
to horizontal. The theoretical predictions and experimental data (each point 
corresponding to 1000 estimates) are compared with the SQL and the Heisenberg 
limit. Both the hybrid and non-adaptive schemes surpass the SQL with increasing 
iV, and are consistent with respective theoretical predictions. The hybrid scheme 
has Aipcst — 0{{\n Ny^^N~^/'^), but the non-adaptive scheme exhibits scaling at the 
Heisenberg limit A0est = 0(A^^^). 



The experimental data and numerical predictions for both schemes are shown in 
figure 131 We vary the number of resources by varying the maximum number of passes 
2^. The limited aperture of our optics (50 mm) lead to a limit of maximum 32 passes 
(maximum K = 5). For this reason, the hybrid scheme had a maximum A^ of 95, and 
the non-adaptive scheme (with Mk = 2 and /i = 3) a maximum A^ of 297. 

The results for the hybrid scheme agree with our theoretical prediction of an 
asymptotic scaling of A(/)est close to 0(A^~^/'^). This shows the interesting result 
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that the SQL can be surpassed by combining two measurements that do not do so 
individually. The data for the non-adaptive scheme are also consistent with numerical 
predictions, clearly surpassing the SQL and exhibiting scaling at the Heisenberg limit 
^0est = 0(A^~^). For the maximum number of resources used, the data demonstrate the 
predicted overhead factor of 1.91 (slightly below the asymptotic factor of 2.03) relative 
to the Heisenberg limit. 

Finally, we make some observations about the role of photon loss in the techniques 
described. With either NOON states or multipass single photons, the interferometric 
scheme is relatively sensitive to loss in the sample [39], with efficiency r/overaii oc r]'^ , 
where t] is the efficiency for a pass of a single photon through the sample. However, the 
multipass scheme is less sensitive to loss at the detectors, where ^^overaii oc r], as opposed 
to r/overaii oc T]^ foY NOON states. To reduce the effect of loss, it may be valuable to 
employ phase-sensitive states that are loss-tolerant, such as those described in [38| 139]. 
It should be possible to configure both adaptive [30] and non-adaptive schemes for use 
with these states. 

6. Conclusions 

When using an n-photon NOON state, or n passes of a single photon, in a two-mode 
interferometer with an unknown phase in one mode, the output probabilities vary 
sinusoidally with n(f) rather than (f). Although this gives phase resolution of 0(l/n) in 
principle, it does not give a phase estimate with this accuracy, because it is sensitive 
only to changes modulo 27i/n. This ambiguity poses a fundamental problem, because 
the accuracy of the prior information about the phase required to resolve the ambiguity 
is of the same order as the phase resolution obtained by the measurement. In this 
paper we have devised two simple yet highly efficient techniques for eliminating the 
phase ambiguity. For the second technique, we obtain an unambiguous phase estimate 
with accuracy scaling as 0(1/A^) in the total number of photon-passes A^. This is the 
Heisenberg limit scaling, in that it is the best possible under the restriction of the 
uncertainty principle for N and 0. Moreover, we have experimentally demonstrated our 
schemes using multiple passes of single photons. 

In contrast to the methods presented in [H |30], the Heisenberg-limited 
scheme introduced here does not require adaptive measurements, but instead 
uses a predetermined sequence of measurements. Reducing or removing adaptive 
measurements should make practical implementation of precision phase estimation 
far easier. The crucial feature which enables 0(1/A^) scaling with nonadaptive 
measurements is the use of larger numbers of repetitions for the measurements with 
fewer passes. Using this, we have analytically proven that 0(l/iV) scaling is obtained. 
In contrast, it is not possible to obtain the Heisenberg limit scaling with nonadaptive 
measurements if the number of repetitions is held constant. 
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Appendix A. Quantifying resources 

In this work we quantify the resources by the total number of applications of the phase 
shift, regardless of whether they are applied in parallel (to different photons) or in 
series (repeatedly to the same photon). This allows the resources to be treated in a 
consistent way, yielding a theoretical limit scaling as 0(1/A^) [3]. One might argue that 
the resources could be quantified by a pair of numbers: the number of photons and the 
number of passes. Then an n-photon NOON state with a single pass and n passes of a 
single photon would be regarded as using different resources. However, this distinction 
cannot be maintained rigorously, because the experiment can be designed such that 
exactly the same quantity is being measured in either case. We now demonstrate this 
fact. 

Consider the case that the NOON state is detected by time-resolving photon 
counters. For a two-mode state of sufficient duration, the timings of the detections 
will typically be distinct. The time that each photon is detected can be used to 
determine the time that it passed through the phase shift. (For the experiment we 
consider superpositions of left- and right-circular polarization, so the photon always 
passes through the phase shift, in contrast to the case where different spatial modes are 
used.) Given a set of arrival-time data {tm} from a NOON-state experiment, one can 
design a single-photon experiment (using a set of mirrors and electro-optic switches), 
such that the single photon passes through the sample at the same times tm- In each case 
the phase shift is sampled in exactly the same way, by photons passing through at times 
tjn- Because the photons are indistinguishable, it makes no difference whether it is n 
photons in a NOON state, or the same photon passing n times. Similar ideas have been 
explored in [U ISl |4l [23l EU |29] . Interestingly, it is not difficult to devise a series of thought 
experiments that interpolate between the NOON state experiment and the multiple-pass 
experiment. For example, consider the case where the photon passes through the sample 
n times, but it is duplicated between passes. This duplication happens in a coherent 
way, so the state a|0) + (3\1) becomes a|00) + /3|11), where indicates left-circularly 
polarized, and 1 indicates right-circularly polarized (see figure lAT!) . The interpretation 
is that the lower photon in the spht is the "same" photon, so the photon passes n times 
through the sample. On the other hand, because the photons are indistinguishable, we 
could equally interpret the upper photon in the split as the "same" photon. Then we 
have n different photons, each passing through the phase shift once. This experiment can 
then be transformed to the NOON state experiment simply by moving all duplications 
of the photons to the beginning. 

The difference between n passes of a single photon, and single passes of n different 
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Figure Al. A thought experiment intermediate between phase measurements with a 
NOON state and multiple passes of a single photon. The horizontal lines are photons, 
and the splittings indicate that the photon has been duplicated coherently. 

photons is simply due to the arbitrary interpretation of which photon is the "same" 
photon. This therefore makes it ambiguous as to how the resources should be quantified 
if the photons and passes are quantified as different resources. Another difficulty 
with that way of quantifying resources is that in a complicated measurement involving 
different photons with different numbers of passes, it is ambiguous as to which number 
of passes should be used. In contrast, the number of passes of single photons through 
the phase shift is always unambiguous, and has a rigorous mathematical foundation [3] . 
We emphasize that, although non-classical states and multiple passes may be regarded 
as equivalent, they present different experimental challenges, with different advantages 
and disadvantages. 
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